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Abstract 

The multiple quantum (MQ) NMR dynamics in the system of equivalent spins with 

■ the dipolar ordered initial state is considered. The high symmetry of the Hamiltonian 

5^ . responsible for the MQ NMR dynamics (the MQ Hamiltonian) is used in order to develop 

Qh| the analytical and numerical methods for an investigation of the MQ NMR dynamics in 

' the systems consisting of hundreds of spins from "the first principles". We obtain the 

\& • dependence of the intensities of the MQ NMR coherences on their orders (profiles of the 

MQ NMR coherences) for the systems of 200 - 600 spins. It is shown that these profiles 

i— i' may be well approximated by the exponential distribution functions. We also compare 

43 ■ the MQ NMR dynamics in the systems of equivalent spins having two different initial 

i , states, namely the dipolar ordered state and the thermal equilibrium state in the strong 

^ ' external magnetic field. 

<D '■ 
43 ■ 
O 

The multiple quantum (MQ) NMR dynamics in solids [T] is extremely usefull for an investigation 
of solid structures and dynamical processes therein, for a counting the number of spins in 
Oh" impurity clusters [2J [3], for a simplification of the standard NMR spectra [I]. Usually the MQ 
NMR experiments deal with samples where the nuclear spin system is initially prepared in 
the thermal equilibrium state in the strong external magnetic field pQ. However, it is possible 
qq ! to carry out the MQ NMR experiments with samples prepared in different initial states [5]. 

In particular, one can prepare a spin system in the dipolar ordered state [6] using either the 
^ adiabatic demagnetization method in the rotating reference frame (RRF) [6j [7] or the two- 
^ ■ pulse Broekaert-Jeener sequence [6j Ej. The MQ NMR dynamics with this initial state in small 
O ■ spin systems has been simulated in refs. [HI [10] . Using the dipolar ordered initial state in the 
MQ NMR experiment one should expect the earlier appearance of multiple spin clusters and 
correlations in comparison with the MQ NMR experiment with the thermal equilibrium initial 
state in the strong external magnetic field. In fact, the analysis of the MQ NMR experiments in 
the six-eight spin systems [HI [10] demonstrates that the six-eight order MQ coherences appear 
d • earlier in the experiment with the dipolar ordered initial state. 

One of the basic problems for the theoretical description of the MQ NMR experiments is 
an exponential growth of the density matrix dimensionality with the increase in the number of 
spins. Therefore the modern numerical methods have been developed for simulation of the MQ 
NMR dynamics, which are based either on Chebyshev polynomial expansion [11] or quantum 
parallelism [12]. However, these methods allow one to study the MQ NMR dynamics in systems 
of no more than several tens of spins. A significant progress in this direction has been achieved 
in simulation of the MQ NMR dynamics in the system of equivalent spins [TBI PT41 IT5] , which 
may be prepared, for instance, filling the closed nanopore with the gas of spin-carring molecules 
(or atoms). The matter is that, due to the special symmetry of the Hamiltonian governing the 
dynamics in such spin system, it becomes possible to study the MQ NMR dynamics in systems 
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of hundreds of spins and even more p31 03] ■ The nature of the above mentioned symmetry 
can be clarified as follows. As far as the characteristic time between two successive collisions 
with the nanopore walls is several orders less that the time of mutual flip-flops of any two 
nuclear spins (which is defined by their dipole-dipole interaction (DDI)) [131 E]> it seems to be 
reasonable to use the averaged DDI, which may be obtained by averaging over the spin positions 
in the nanopore [13]. This means that the constant of the averaged DDI remains the same for 
any pair of spins in the nanopore [13], so that the nuclear spins become equivalent. For this 
reason, the Hamiltonian of the nuclear spin DDI in the nanopore commutes with the operator 
of the square of the total spin angular momentum I 2 [151 HH]- Thus, it becomes possible to 
use the basis of the common eigenfunctions for the operator of the square of the total spin 
momentum J 2 and of its projection I z on the direction of the external magnetic field instead 
of the standard multiplicative basis of the eigenfunctions of I z , which yields an exponential 
growth of the Hilbert space dimensionality with the increase in the number of spins p~5[ [T6] . 
As a result, we simplify calculations which allows us to succeed in both the investigation of 
the MQ NMR dynamics in systems of 200 - 600 spin- 1/2 particles [15] and in the study of the 
dependence of the coherence relaxation time on the MQ NMR coherence order and the number 
of spins [17]. Emphasize that the nuclear spin system with the thermodynamic equilibrium 
initial state in the strong external magnetic field is used in refs.[15l [161 fl?] - 

The MQ NMR dynamics in the large system of equivalent spins with the dipolar ordered 
initial state is studied in the present paper. The theory of MQ NMR dynamics of equivalent 
spins with this initial state is given in Secf21 The dependence of the MQ NMR coherence 
intensities on the coherence orders (the profiles of MQ NMR coherence intensities) for systems 
of 200-600 spins is represented in Sec]31 The MQ NMR dynamics in systems with the dipolar 
ordered initial state is compared with the dynamics in systems with the thermal equilibrium 
initial state in the strong external magnetic field in SecJU The basic results are collected in 
SecEl 

2 The MQ NMR coherence intensities in systems of 
equivalent spins prepared in the dipolar ordered initial 
state 

We consider the system of equivalent spin- 1/2 particles with the dipole-dipole interaction (DDI) 
in the strong external magnetic field. The secular part of the DDI Hamiltonian [6] reads: 

Hdz = Djk(2Ij z Ik z — Ijxhx — Ijyhy), (1) 

where = 5-x"(l — 3 cos 2 #•,■&) is the constant of DDI, 7 is the gyromagnetic ratio, is the 
distance between the j'th and fcth spins, 6jk is the angle between the internuclear vector 
and the external magnetic field B , and Ij a (a = x, y, z) is the jth spin projection operator on 
the axis a. Using either the adiabatic demagnetization in the rotating reference frame [HI [7] or 
the Broekaert-Jeener two-pulse sequence [HI [8] one can prepare the spin system in the dipolar 
ordered initial state with the following density matrix: 

p(0) = i exp(-/3if & ) « - 0H dz ), (2) 
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Figure 1: The basic scheme of the four period MQ NMR experiment. 



where — is the inverse spin temperature [k is the Boltsman constant, T is the temperature), 
Z = Tr exp(— (5Hd z ) is the partition function and N is the number of spins. 

If the system under consideration consists of the spin earring molecules (atoms) in the 
closed nanopore, then DDI is averaged (incompletely) by the fast molecular diffusion so that 
the constants of DDI for any spin pair equal each other, i.e. Dj k = D [131 E]. As a result, the 
Hamiltonian ([1]) may be written as follows 



H, 



d~ 



D . „ 



(3) 



where I a = J2i=i hai ( a = x i Vi z )> I 2 = + ^ is t ne square of the spin angular momentum. 

It is important to justify that the high temperature approximation (j2j) is applicable to the 
system of equivalent spins. The simple analisis [H] demonstrates that approximation (j2J) for 
the system with the Hamiltonian ([3]) is valid if 



(3DN < 1. 



(4) 



The systems considered in our work have N = 200 — 600, so that the condition (J4j) is satisfied. 

The MQ NMR experiment consists of four basic periods shown in FigJTJ the preparation 
period r, the evolution period t, the mixing period r and the detection. The MQ NMR 
coherences are generated on the preparation period due to the irradiation of the sample by the 
multiple eight-pulse sequence of the resonance pulses P . Let us express MQ NMR coherences 
in terms of the density matrix of the preparation period. For this purpose note that the 
averaged Hamiltonian Hmq (nonsecular two-spin/two-quantum Hamiltonian [1]) describing the 
MQ dynamics in the system of equivalent spins during the preparation period in the rotating 
reference frame may be written as follows [151 EE] : 



H 



D 



MQ 



{(/+) 2 + (/-f}, 



(5) 



where I + and I~ are the raising and lowering operators (I^ = I x ±il y ). In order to investigate 
the MQ NMR dynamics, one has to find the density matrix p(r) for the spin system solving 
the Liouville equation [6]: 



.dp 

''Tr 



[ h mq,p(t)}, 



(6) 



with the initial condition p(0) = Ha z . This initial condition is obtained from eqs. (j2J) and 
([3]) by dropping both the unit operator and the factor (— /3/2 ), which are not significant for 
the MQ NMR dynamics. Taking into account the Hamiltonians on the different periods of the 
MQ NMR experiment (these Hamiltonians are shown in Fig{T]) one can write the expression 
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for the dipolar energy (Hd Z )(T,t) after the mixing period of the MQ NMR experiment (FigJTJ) 
as follows: 

y(r)e^^ + (r)e^II(r)4} m 

{H dz ){T,t) = Tr | | = (<> 

Tr{e- tA ^-p(r)e iA ^p(r)} 

where p(r) = U{r)Hd z U + (t) is the solution to Eq.(jS]) and U(r) = exp(— iHmqt). It is 
convenient to represent the solution to Eq.fjSJ) as the series [19] 



p(r) = Y,Pk(r), (8) 

k 

where p^ir) obeys the relationship [I z , Pfc(r)] = kp k (r) which can be considered as the definition 
of Pfe(r) (in other words, Pfc(r) collects those entries of the density matrix p(r) which are 
responsible for the fcth order MQ NMR coherence). Then eq.([7j) reads 

(H dz )(r) = ^ e -^ Tr{p.(r)p_ fc (r)} = ^-^^ ( 9 ) 

k LI \ n dzi k 

where the fcth order MQ NMR coherence intensity J k is defined as follows [19] : 

= Tr{ Pfc (r)p ,(r)} _ 

This formula will be used in the numerical analysis of the MQ NMR coherences in Secj3j 

It is easy to write the explicit expression for Tr{H% z }. First of all, using eq. ([3]) we may 
write 

Tr{^} = ^Tr(3/ Z 2 - J 2 ) 2 = (11) 

^Tr{9/ 4 - 3/ 2 / 2 + J 2 (J 2 - 3/ 2 )} = 

3D 2 3D 2 

^Tr{3/ 4 - P Z P} = ^Tr{/ 2 4 - J 2 / 2 }. 

Here we take into account that 

Tr{/ 2 (/ 2 - 3/ 2 )} = 0, Tr{J 2 / 2 } = Tr{/J/ 2 }. (12) 

Next, it is simple to obtain the explicit expressions for Tr{/ 4 } and Tr{J 2 / 2 }: 

Tr{/ 4 } = 2 N ~ 4 N(3N - 2), Tr{/ 2 / 2 } = 2 N ~ A N 2 . (13) 

Finally, eqs. flimiUl) yield the following result: 

Tr{Hj z } = 3N(N - 1)2 N ~ 4 D 2 . (14) 

Taking into account the structure of the MQ Hamiltonian Hmq, it may be readily shown that 
only the even order MQ NMR coherences appear in our numerical experiment and the coherence 
order may not exceed the number of spins N [TJ [T9~] . 
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It is obvious that the MQ Hamiltonian fl5]) commutes with the square of the total spin 
momentum operator I 2 [151 HSJ 02] • As f ar as t-^ 2 ' = 0' ^ * s possible to use the basis of the 
common eigenf unctions of operators I 2 and I z for the description of the MQ NMR dynamics. 
Namely this fact allows one to avoid the problem of the exponential growth of the matrix 
dimensionality with the increase in the number of spins, which appears in the traditional 
multiplicative basis [T] of the eigenfunctions of the operator I z . In the new basis, the MQ 
Hamiltonian Hmq consists of blocks H M q corresponding to the different values of the total 

spin momentum S (S 2 = S(S + 1), S = y, y — 1, y — 2, ■ • • , y — [y] , [i] is an integer part of 
i): 

- JY_ JY_ i _ JV_[JV] 

^ A/Q = diag{^ Q ,^- ,...,ff£ g laJ }. (15) 

As far as the Hamiltonian //^ (J3J) is diagonal in the basis of the common eigenfunctions of 
the operators I 2 and I z , then the density matrix p(0) may be also splitted into the diagonal 
blocks p s (0) with S — y , y — 1, y — 2, ■ • • , y — [—1 . Consequently, the matrix p(r) has also 
the diagonal block structure with blocks p s {r). Thus, the problem becomes separated into the 
set of independent problems for each (2S + 1) x (2S + l)-dimensional block p s (r) which is a 
solution to the Liouville equation ([6]) with the Hamiltonian H M q- Of course, expansion (|8]) 
may be applied to each block p s . The contribution J^s from the block p s to the intensity Jk 
of the /cth order coherence is defined by the obvious formula j!5j : 



T ( , Tr{pf(r)pg fc (r)} 

JkAr) = ^ , (16) 

where pf (r) is the contribution from the matrix p s to the fcth order coherence. One has to take 
into account that each block p s {r) is degenerated with the multiplicity n N (S) [201 115] : 

= ( f + s + i)!(f- s) r T' (17) 

As a result, the observable intensities Jfc(r) (— < k < N) are following [T5l fT6] : 

^(T)^%(5)J fe , s (r). (18) 

5 

Remember that the dimensionality of each block H M q of the MQ NMR Hamiltonian Hmq is 
(2S' + 1). Taking into account the block degeneration we obtain the correct value for the matrix 
dimensionality of both the Hamiltonina Hmq and the density matrix p [15j : 

J2n N (S)(2S+l)=2 N , (19) 

s 

which is valid for the system of iV interacting spin-1/2 particles. 

The numerical algorithms describing the MQ NMR dynamics in the systems of equivalent 
spins with the thermal equilibrium initial state in the strong external magnetic field have been 
developed in [151 EH 02] • With minor corrections, these algorithms may be used for the simu- 
lation of the dynamics of the MQ NMR coherences in the spin system with the dipolar ordered 
initial state. In particular, the integral of motion related with the MQ NMR Hamiltonian in- 
variance with respect to the rotation over the angle n around 2;-axis [21] is also present. Thus, 
for the odd N, it is enough to solve the problem for the MQ NMR Hamiltonians with two 
times lower matrix dimensionality and then double the resulting intensities [21]. We use the 
spin systems with odd N in all numerical calculations in the next section. 
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3 The numerical analysis of the MQ NMR profiles 



Using the method developed in the previous section we investigate the profiles of MQ NMR 
coherences. As far as all spins are "nearest neighbors" in the system of equivalent spins, N 
spin cluster appears already after the time interval r ~ 1/D. However, some reorganization 
of this cluster is required for the MQ coherence formation [15]. The analysis of the MQ NMR 
coherence dynamics demonstrates that the quasistationary profile of MQ NMR coherences is 
created during f ~ 2 (f = Dt is the dimensionless time hereafter) and remains fast oscillating 
for r > 2. Because of these oscillations it is convenient to use the averaged intensities Jk 
[T5] instead of intensities Jk itselves. We estimate the dimensionless averaging time interval 
as T ~ 27r/|A™ 2 n | 7.255, where A™ 2 n = y/3/2 is the minimal eigenvalue of the Hamiltonian 
[T5l IT6] . For our convenience, we take T = 8 so that 

- 1 [ 2+T 

Jk = f J Jk{f)df. (20) 

It is observed that Jk do not significantly vary with increase in T, so that the definition of 
the averaged intensities given by eq. (T20|) is valid. Although the dynamics of all coherences 
has been found in the numerical simulations, the intensities of the high order coherences are 
negligible, so that we represent the intensities of the MQ NMR coherences up to the 50th order 
in all figures below. The profiles of MQ NMR coherence intensities for systems of 201, 401 and 
601 spins with dipolar ordered initial state are shown in Figf2J These profiles are similar to 
those which have been found for the systems with the thermal equilibrium initial state in the 
strong external magnetic field [151 US]- Similar to refs.|15[ IT6] . the averaged intensities of MQ 
NMR coherences are separated into two families: 

ri = {J 4fc _ 2 , fc = 0,±l,±2,...}, (21) 

r 2 = {J 4fc , fc = ±i,±2,...}, 

with the zero order coherence intensity Jo does not corresponds to any of these families. Each 
family may be approximated by the smooth distribution function as follows: 

4 

A 1 (l + ^(-l)'a l4 (2|A;|r)e- 2 " 1 l fc l, k = ±1, ±3, . . . 

l f , (22) 

A * (l + J](-l) J a 2j (2|A;|) J )e- 2Q2 l fc ', k = ±2, ±4, . . . 



J-lk ~ < 



i=l 



where parameters A4, and a\ for spin systems with N = 201, 401 and 601 are collected in 
Table 1. The algorithm defining the approximation parameters is similar to that suggested in 
ref.[T2]. 

Similar to the profiles of the MQ NMR coherence intensities obtained for the systems of 
equivalent spin-1/2 particles with the thermal equilibrium initial state in the strong external 
magnetic field [15] , profiles for the systems with the dipolar ordered initial state seemed out to 
be exponential. This conclusion agrees with results obtained during an elaboration of numerous 
MQ NMR spectra [22] and contradicts the phenomenological theory [1] predicting the Gauss 
profiles of the MQ NMR coherence intensities. 
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N 


M 


an 


a 12 


a 13 


ai4 


«i 


201 
401 
601 


8.498 x 10~ 2 
4.546 x 10~ 2 
3.546 x lO- 2 


3.133 x lO" 1 
1.694 x 10" 1 
1.456 x 10" 1 


8.025 x 10~ 2 
3.357 x 10~ 2 
2.471 x lO" 2 


6.109 x lO" 3 
1.649 x lO" 3 
1.028 x 10" 3 


2.739 x lO" 4 
5.739 x 10- 5 
2.855 x 10" 5 


2.437 x 10- 1 
1.638 x lO" 1 
1.342 x 10" 1 


N 


z 


aoi 


CJ09 


Oo-i 
w zo 


a 94 


dor. 


201 
401 
601 


2.217 
1.411 
1.210 


4.029 x 10" 1 
3.587 x 10" 1 
3.441 x lO" 1 


7.522 x lO" 2 
5.648 x lO" 2 
5.123 x lO" 2 


7.638 x 10" 3 
4.616 x lO" 3 
3.924 x lO" 3 


4.624 x 10" 4 
2.116 x lO" 4 
1.664 x lO" 4 


1.562 x 10" 5 
5.136 x lO" 6 
3.686 x 10" 6 


N 


0,26 






201 
401 
601 


2.510 x 10~ 7 
5.533 x 10~ 8 
3.512 x 10" 8 


3.039 x lO" 1 
2.183 x 10" 1 
1.940 x 10" 1 



Table 1: The parameters A iy a>ij and «j of the distribution function given by eq. (|22~l) for 
N = 201, 401 and 601 



Liouville equation fl6]) and formulas for the MQ NMR coherence intensities ( fit)]) yield the 
conservation law of the sum of the coherence intensities [23] : 

£J*(t) = 1. (23) 

k 

This law together with the approximating formula fl22|) allows one to find a good approximation 
to the zero order coherence intensity. We compare the calculated values of Jo (see eq.f JTO]) ) with 
the values Jq PW = 1 — 2 Y^kLi f° un d using the above conservation law and distribution 
function ([22]): 

( 2.519 x HT\ iV = 201 ( 2.501 x 10 _1 , ^ = 201 

J = < 2.417 x 10 _1 , N = 401 , J appr = <^ 2.381 x 10 _1 , = 401 . (24) 
{ 2.361 x 10- 1 , A^ = 601 [ 2.387 xlO- 1 , = 601 

Some discrepancy between Jo and Jq PW appears because we take into account only coherences 
up to the 50th order in constructing the distribution function ( |22|) . while contributions from 
the higher order coherences are missed. 

4 Comparison of the MQ NMR dynamics in the spin 
systems with two different initial states 

The preparation of the system in the dipolar ordered initial state means that the two-spin 
correlations appear already at the initial time, unlike the standard MQ NMR experiment, 
where the thermal equilibrium initial state is defined by the one-spin Zeemann interaction with 
the external magnetic field (9j [10]. This statement is justified by Figj3j where the formation 
times ff(n) of different order coherences for both initial states are represented. Here the term 
formation time ff(n) of the nth coherence means the time moment when the nth coherence 
intensity J n (f) reaches the value J n for the first time. We see that MQ NMR coherences in 
the system with the dipolar ordered initial state (the lower solid line) appear much earlier. 
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Coherence order, n 



Figure 2: Coherence intensities profiles for spin systems with N = 201,401,601. The inset 
demonstrates that zero order coherence does not belong neither to Fx nor to IV Only intensities 
of positive order coherences are presented. 



i -, 




Coherence number, n 

Figure 3: Coherence formation time f/ versus the coherence number for the dynamics of equiv- 
alent spins with the dipolar ordered and the thermal equilibrium in the strong external field 
initial states. 
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Coherence order, n 

Figure 4: The MQ NMR profiles for the systems of equivalent spins with the dipolar ordered 
and thermal equilibrium in the strong external magnetic field initial states, N = 201; the inset 
demonstrates that Jo and J4 are essentially bigger in the case of the dipolar ordered initial 
state. 

This result agrees with that obtained in [9] for the MQ NMR in systems with small number of 
spin-1/2 particles. 

The profiles of the MQ NMR coherence intensities for the system of N = 201 spins with 
both the dipolar ordered and the thermal equilibrium in the strong magnetic field initial states 
are compared in FigJH This figure demonstrates that the discrepancy between two families of 
MQ NMR coherences I\ and T 2 (see Eq. (l2ip ) is bigger for n < 10 and smaller for n > 10 in the 
case of the dipolar ordered initial state. The inset shows that J and J 4 are essentially bigger 
in the case of the dipolar ordered initial state. 

Thus, similar to the usual NMR experiments in solids [B], the MQ NMR experiment in the 
systems of equivalent spins with the dipolar ordered initial state may be usefull to supplement 
the MQ NMR experiment with the thermal equilibrium initial state in the strong external 
magnetic field. 

5 Conclusions 

We have studied the MQ NMR dynamics in the systems of equivalent spin-1/2 particles with 
the dipolar ordered initial state. For this purpose we modify the method developed in ref.[T5] 
for the system of equivalent spin-1/2 particles with the thermal equilibrium initial state in the 
strong external magnetic field. Similar to ref.[15], the high symmetry of such systems allows 
one to investigate the dynamics in large spin systems containing hundreds of interacting spins. 
We obtain dependence of the MQ NMR coherence intensities on their order (the profiles of the 
MQ NMR coherence intensities) in systems of 200-600 spins and demonstrate that these profiles 
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may be well approximated by the exponential distribution functions. As far as the analogues 
result has been obtained in refs.|15[ [16] for the systems with the thermal equilibrium initial 
state in the strong external magnetic field, we may suppose that the exponential profiles of 
MQ NMR coherence intensities are fundamental fact in MQ NMR dynamics. The theoretical 
results obtained in refs. [24"| [25] confirm this conclusion. 

We demonstrate that the MQ NMR coherences appear faster in the spin systems with the 
dipolar ordered initial state. The MQ NMR experiments with the dipolar ordered initial states 
expand possibilities of the MQ NMR spectroscopy in study of the structures of solids and the 
dynamical processes therein. 

All numerical simulations have been performed using the resources of the Joint Supercom- 
puter Center (JSCC) of the Russian Academy of Sciences. The work was supported by the 
Program of the Presidium of Russian Academy of Sciences No. 21 " Foundations of fundamental 
investigations of nanotechnologies and nanomaterials" . 
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